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2. 

1,  Introdaction. 

The  proTjlem  of  propogation  of  radio  waves  in  a  homogeneous  at- 
mosphere has  "been,  studied  theoretically  "by  several  iavastigators,  "begin- 
ning atout  1909  with  the  v;ork  of  SoEaierfeld,  A  "brief  account  of  this  work 
will  he  found  in  Stratton's  Slectroaegnetic  theory,  p. 573  ff.  In  this 
work  the  effect  of  reflection  from  the  earth  and  from  the  ionosphere  has 
heen  treated.  This  problem  can  he  said  to  "be  largely  solved,  at  least 
for  long  wave-lengths. 

In  the  case  of  an  inhomogeaeous  non-standard  ataiosphere  the 
prohlera  of  propagation  offers  great  difficulties.  Attempts  have  "been 
made  to  approximate  actual  atmospheres,  and- the  resulting  pro'blems  have 
"been  solved,  also  approximately.  One  such  method  -  that  carried  out  hy 
v;,  H.  Furry  -  approximates  some  possi"ble  atmospheric  condition  hy  two 
layers,  each  v;ith  a  constant  refractive  index.  The  method  is  theoret- 
ically suitable  for  any  number  of  layers.  In  this  study  the  layers  are 
plane,  and  in  order  to  obtain  results  applicable  to  a  spherical  earth 
surrounded  by  layers  in  the  form  of  spherical  shells,  an  earth  flatten- 
ing factor  is  used.  This  is,  of  course,  a  further  approximation.  Jilven 
the  approximate  solutions,  for  limited  cases,  have  led  to  immense  amounts 
of  calculation  in  which  progress  has  been  slow. 

The  present  paper  represents  an  attempt  to  deal  with  the  prob- 
lem of  propogation  in  a  special  case  of  an  ichomogeneous  atmosphere  so 
as  to  obtain  an  exact  solution.  The  problem  is  the  following:  We  con- 
sider a  spherical  earth,  v/hich  is  finitely  conducting  and  within  which 
the  index  of  refraction  is  constant,  surrounded  by  a  layer  of  finite 
depth  in  which  the  index  of  refraction  varies  inversely  as  the  square  of 
the  distance  from  the  earth's  center.  Outside  this  layer  is  an  infinite 
layer  in  which  the  refractive  index  is  constant.  A  source  of  radiation 
is  situated  in  the  finite  layer,  ^e  are  to  find  the  electromagnetic 
field  at  any  point. 

We  have  restricted  the  layer  with  variable  index  of  refraction 
to  finite  depth  since  it  is  clear  that,  whilo  such  a  variation  of  the 
refraction  index  Is  possible  in  a  finite  layer,  if  the  layer  were  ex- 
tended infinitely  far  we  would  have  the  refractive  index  tending  to  0, 
which  is  physically  untrue. 


3. 


Ihe  method  of  attaclc  is  the  follov/ing.  The  wave  aquation  In 
spherical  coordinatoe  satisfied  hy  the  Hertzian  potential  function  in  the 
case  of  a  constant  index  of  refraction  is 

V^  u  ->  k^u  =  0.  (A) 

A  coiqplete  set  of  solutions  of  this  equation  is  given  in  teros  of  Bessel 
and  Hankel  functions  and  the  Legondre  "Doiynoiiiials ,   in  the  case  of  syncao- 

3.2 

try  about  an  axis.  If  ue  let  r  =  e  where  a  =  the  radius  of  the 

r 

earth,  (this  is  the  inversion  in  the  sphere  representing  the  earth's  sur- 
face, with  the  earth's  center  as  the  center  of  inversion)  equation  (A) 
tocoiaes  «  p 

V^U  <■  ^-g-U-  0  ,  (B) 

.  ,  r 

and  it  can  he  shotm  that  the  fu-action 

u(r.a,  cp)  =  i  u<-|—  ,  e,  95) 

satisfies  this  modified  wave  equation.  Hence,  from  the  complete  set  of 
solutions  of  (a)  we  are  ahle  immediately  to  vnrite  down  a  complete  set  of 
solutions  of  (B),  Since  (B)  is  the  wave  equation  satisfied  hy  the  Hertzian 
potential  function  in  the  case  of  a  refractive  indez  which  varies  in- 
versely as  the  square  of  r,  we  are  in  a  position  to  solve  the  problem 
described  above,  provided  the  boundary  conditions  can  be  met. 

The  boxindary  conditions  are  derived  from  Idaswell's  equations. 
In  the  case  we  are  considering  the  quantity  6  ,  usually  called  the 
dielectric  "constant",  is  not  constant  but  varies  inversely  as  the  fourth 
■Dower  of  r.  It  is  found  that  the  boimdary  conditions  in  this  case  are 


essentially  the  seme  as  for  a  constant  £  .  These  conditions  are  that  £u 

and 

layer. 


and   — ^ — ~    must  be  contimioas  across  the  bounding  surfaces  of  the 
&  r 


The  solution  of  the  problem  is  straightforward  and  only  minor 
approximations  have  been  used.  The  result  is  in  the  form  of  am  infinite 
series  of  Bessel  and  Hankel  functions,  with  coefficients  found  from  four 
simultaneous  linear  eqtiations.  If  the  series  converged  rapidly  the  cal- 
culation of  the  necessarj'  coefficients  would  not  be  an  insurmountable 
obstacle.  But,  unfortunately,  the  series  converges  very  slowly  and  the 


K 


work  necessary  to  obtain  useful  results  would  "be  prohibitive.   It  must  "be 
remembered  that  results  calling  for  a  formidable  amount  of  calculation 
have  been  encountered  before,  notably  in  the  case  of  Watson's  solutions 
of  the  propagation  problem  in  a  homogeneous  atmosphere.  The  work  of 
Eckersley*  and  that  of  van  der  Pol  and  Bremmer**  was  aimed  at  simplifica- 
tions of  these  results  so  as  to  make  the  calculations  reasonably  simple. 
It  is  conceivably  that  the  results  obtained  in  this  paper  can  also  be 
brought  within  the  range  of  practical  coniputation.  At  present  it  does 
not  seem  that  the  case  covered  has  sufficient  application  to  warrant  the 
attempt,  at  least  not  until  other  avenues  have  been  explored. 

This  work  has  some  relationship  to  that  of  W.  H.  I'uTTy   men- 
tioned above,  as  may  be  seen  from  the  following  consideration.  Parry 
makes  use  of  the  modified  index  of  refraction,  iJ,  defined  as  follows: 

„  _  (index  of  refraction)  x  (distance  to  earth's  center) 
same  product  taken  at  a  special  reference  point 

In  the  case  of  two  layers,  each  v/ith  a  constant  index  of  refraction,  the 
modified  index  has  the  f oim 


H  =  m^r. 


r  <  r-. 


r  >  r-, 


where  r     is  the  value  of  r  at  the  interface  between  the  layers.     Ihis 
leads   to   the  so-called  bilinear  M-curve  made  lip  of  tvfo  segments  with 
different  slopes  as  shown  in  Tig.   1.     With  two  layers,   in  the  first  of 


/I. 


A. 


A/ 


/=/c^   1 


// 


/^/<^   -? 


*  Sckersley  and  Killington,  Royal  Soc.  Phil.  Trans.,  v.  237,  1939,  P  273; 
Phil.  Ma£..v.27,  1939.  P  517. 

**  van  der  Pol  end  Bremmer,  Phil.  I^.ag. .  v.  2U,  1937,  P  1^1,  P  S25:  v.  25, 
193s.  p  SI7. 
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which  the  index  of  refraction  varies  inversely  as  the  square  of  r  and  in 
the  second  of  which  it  is  constant,  we  have  for  the  modified  index 

K  =  ^\  r  <  r^. 

N  =  Digr,  r  >  r^. 

The  M-curve  in  this  case  is  that  shown  in  Fig.  2. 

It  should  he  emphasized  that  we  are  dealing  with  a  spherical 
earth  from  the  "beginning  and  no  earth  flattening  device  is  necessary. 
This  study  raay  he  considered  to  he  a  first  step  in  an  attempt  to  solve  the 
propogation  problem  for  inhonogeneous  atmospheres  closely  approximating 
true  conditions.  It  seems  hest  to  follow  up  this  work  hy  the  use  of  other 
expressions  for  the  index  of  refraction  that  \d.ll  leave  the  wave  equation 
in  some  manageable  form  and  will  leave  the  houndarj'  conditions  essentially 
the  seme  as  those  used  here.  An  attractive  possibility,  which  will  be 

investigated,  is  that  in  which  £-  =  i  +  _£ —   where  r  is  the  distance 

r^ 
from  the  center  of  the  earth.  Under  this  assumption  it  v;ill  not  be  nec- 
essary to  consider  the  atmosphere  to  be  in  layers;   £.  can  vary  continu.- 
ously  from  r  =  radius  of  the  earth  to  r  =  infinity,  and  as  r  tends  to  in- 
finity S  tends  to  1  as  physical  considerations  show  that  it  should.  In 
this  case  the  only  boundary  surface  oa  which  special  conditions  will  have 
to  be  met  will  be  the  earth's  surface,  ■* 

2,  Mathematical  Formulation  of  the  Problem. 

Since  we  are  interested  in  propagation  in  a  layer  in  which  & 
is  variable,  we  must  examine  the  wave  equation  to  be  satisfied  in  this 
layer  and  also  the  botmdary  conditions  that  must  be  met. 

Maxvrell'a  equations  are 

VXE=--||.      (1)  '     ^XH  =  -|f.        (2) 
V  B  =  0,        (3)       V'  D  =  0,         (U) 


in  which  D  =  £E  and  B  =/^H,  Ve  assume y6<  to  be  constant  but  suppose 
that  £  is  variable  from  point  to  point.  Me   shall  introduce  a  Hertzian 
vector,  TT  ,  in,  the  usual  way. 
Let 

H  =  V  )^-^-^SP-  •  (5) 


6. 


This   expression  autcciatically  satisfies   (3).      Su'bstituting  in  (l), 

2> 


VXS<-/<V)^-^-i^:iiL  =0, 


at2 


SMs  eqiiation  fill  "be  satisfied  if  we  put  the  quantity  in  parentheses 
equal  to  the  gradient  of  a  scalar  function.  Let 

S+/£-^~E-  =^4^  ,  (6) 

uhere  vjy  is  yet  to  'be  detenained. 

Su'ostituting  froa  (5)  and  (6)  into  (2),  uo  havQ 

Because  of  the  nature  of  the  propagation  proTjlem  in  which  W3 
are  interested,  we  adopt  a  cpherical  coordinate  system  and  assume  that  S 
depends  upon  r  only,  For  later  ptirposes  it  will  "be  sufficient  to  assuoie 
that  It  has  only  an  r-coapojient ,  "[!_«  Then 

y*   r  sin©  ^r^       \'       F    ae J'       ^^^ 

and 


i.   • 


1  a^^^^r)      1     d(tirr)\        (g) 


r  d^T    bo        *  r  sia  t)  ar  d  Q? 

Sine 


)■ 


V  M^  Vdr  '  r  "S^'  r  sin  9   a  (p  y* 

vo   find,  using  the  8-  and  (p-  components  as  deterained  fron  (7) 
and  (9), 

r     h  O  r   a  r  od    » 


(10) 


7. 


and  ■  ; 

r  sin  e  &(p  r  sin  ©   a  r  3  0      *  ^^ 

Equations  (10)  and  (11)  are  satisfied  if  we  talce 

Putting  tbis  in  the  r-co2iponent  as  found  from  (7)  and  (9), 
we  have 

or 

We  now  assume  that  If  depends  on  the  time  throu^  a  factor  of 

-icot  ^ 

or 

We  nay  then  write. 


-icot 
the  form  e     only,  where  cois  2n  times  the  frequency  of  the  radio  waves. 


Putting  this  in  (I3)  and  dividing  hy  r,  we  get  the  modified  wave  equation: 


1  ^^(ra)^^aUu)^^,£.£l(eO%     1    ^  U   9  i^A 

^-2-V7|^*/^2^-=°» 
r  sin  Q   <^(r 


or 


"  /r\2 


V^u-H^-^  =u.  ^^i^^^^-^/e^\  =  o.  (Ill) 

tr  d  r  9J^ 


In  the  case  in  which  we  are  interested  u  =  A"^  -,  ■■  and  £=  ^ —  • 


s. 


In  these  ezprossions 

r  =  distance  from  any  point  to  the  center  of  the  earth, 

a  =  radius  of  the  earth, 

k  =  2  TT 

X"  ' 

E  =  distance  from  the  transmitter  to  any  point  in  the  mediiim 


~  /r^  *  1)^  -  2rh  cos  6  , 

xih.ezz'b   is  the  distance  from  the  trsmsmitter  to  the  center  of  the  earth 
and  9  is  the  angle  "between  h  and  r. 

With  these  values,  keeping  in  mind  that  £  differs  hut  little 

d  R 

from  1,    that  R  and     ^       are  small  in  coisparison  with  r,  we  find  that 

'  11 

the   second  and  third  terms  of  equation  (lU)  are  roughly     —  and  — ~ 

tines  the  magnitude  of  the  first  and  last  terms,   respectively.     Ve  see, 
therefore,   that  we  laay  assume  that  u  must  satisfy  the  wave  equation  in 
spherical  coordinates , 

y2  U     ■^/<eo\  =   0.  (15) 

From  the  knowledge  of  u  we  can  find  E  and  H  from  (5)  and  (6), 

-iCJt 

Ve  have,  suppressing  the  factor  e     , 

E^-/<etJ-ru*  ~^|^--g-j      .  Eq  ,-^-g_-__j. 

E  =     1      a^(&ru)  ,,^v 

^<P       f  r  sin  0   dr  dc/*   *  ^  " 

If  we  assume  asial  symmetry,  derivatives  with  respect  to  Cf  vanish.  On 

the  "boundaries,  therefore,  ve  mtist  satisfy  the  conditions  that  Em   and  Eg 

are  continuous.  Since  r  is  constant  on  the  "bounding  surfaces,  we  may 
Integrate  E(p   and  Bq  with  respect  to  6  .  Hence  we  find  that  £.  u  and 

1  _g  V.  fc  ru ;   jnngt  ^.q   continuous  on  the  "boundaries.  The  second  condition 
£   ^  r 

requires  that   -— ^ "^    -|^  ru  "be  continuous,  end,  just  as  in  the  case 

of  oquation  (lU),  ve  can  show  that  -~  ru  may  "be  neglected  in  comparison 

o 


9. 


yith.     — i£iii.  ,     \is  can  tiierefora  make  use  of  the  following  1)0-00.(13X7  condl- 
tions,   naaely,    that     &  u  axid    -^rrr — ^^st  he  continaous, 

3,   Inversion  of  the  Wave  Equation. 

Let  u(  P  ,&,    ^)  satisfy  the  uave  equation, 

y2^       *       ^^u  =  0.  (IS) 

It  can  he  shown  "by   direct  sobstitutioa  that 

U(r.  Q,fP)   =  ^  u(  ~-  ,  O.C?)  .  (19) 

Edtisfies  the  zaodifisd  vsstq  epilation, 

2  2 

V^  *  ^  \  u  °  0,  (20) 

^2 

lato  which  (IS)  is  transformed  hy  puttir^  /^  =  -p-  • 
If  we  pat 

/S~  '^n  +  'It-^   '  (21) 

(1) 

in.  which  E    is  the  Hankel  function  of  the  first  kind  and  J  is  tha 

Bessel  function  of  the  first  kindj  the  functions  "2^  (kP)?  (coa  ©)  and 
\/  (k  P  )  P  (coa  0)  fona  a  complete  set  of  solutions  of  equation  (IS) 

in  the  case  of  axial  syisaetsy.  In  these  functions,  the  quantities 

?  (co3  0)  are  the  Legendre  polynoiaials  with  argument  =  cos  Q.  It  folloirs 

that  the  functions 

fora  a  cosrplete  set  of  solutions  of  equation  (20),  in  the  case  of  asial 

sycnetry, 

U,  romal  Solution  of  the  Proposed  Propagation  Prohlem. 

V/3  now  consider  the  propagation  prohlea  stated  in  the  intro- 
duction. The  earth  is  assuned  to  he  a  homogeneous  inperfectly  conducting 


¥n'-)=./@:  J  .„t) 
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spiiere  of  radius  a.   In  a  layer,  I,  'bounded  Toy  the  earth  and  an  outer 
spherical  surface  of  radius  d,  the  refractive  index  is  assumed  to  vary 
inversely  as  the  square  of  r,  the  distance  to  the  earth's  center.  Oat- 
side  this  layer  the  refractive  indez  is  a  constant  equal  to  that  for  free 
space.   V/e  assume  that  a  transmitter,  hehaving  like  a  Hertzian  dipole  of 
angular  velocity  to,  is  situated  in  I,  at  distance  "b  from  the  center  of 
the  earth.  We  shall  evaluate  the  field  at  any  point  P  with  spherical 
coordinates  (r,  Q,  cp). 

In  order  to  ©"btain  the  solution  of  this  prohleia, which  is  associ- 
ated with  equation  (20),  we  first  consider  the  propagation  prohleai  associ- 
ated v/ith  equation  (IS).  This  is  the  case  in  vAiich  the  refractive  indes  is 
constant  in  the  layer  containing  the  transmitter.  The  field  at  any  point 
in  this  layer  will  "be  the  resultant  of  the  field  of  the  transmitter  and 
the  field  produced  "by  the  earth  and  the  outer  region.  \/e  therefore  put 

•u,  =  u.  +  u  ,  in  which  u  is  the  Hertzian  potential  for  the  transmitter, 
t    a  t 

and  u  is  the  Hertzian  potential  for  the  distur"bed  field.  Letting 
^  ikR 

u,  =   f ,  ^  .  where  R  represents  the  distance  from  the  transmitter, 
t     ikR  ' 

we  have  * 

ikS    CO 
u^  =  V^  =  Y"  (2n  *  1)  r  (khAj/^ (kr)  P„(cos  6),  r  <  h,      (23a) 

\ 

00 

=  XI  (2n*  l)^Jkr)-vy^(kh)  Pjcos  O),        r  >  Ij,      (231j) 

where  '^  and  'vl/'  have  the  meanings  assigned  to  them  in  (21)  and  P  are 
the  Legendre  polynomials.  Me  may  also  write 

u  =  T"  (2n-»-  1)  r^  'a/_(kr)  +  c_  r„(kr)l  P_(cos  6),     (2U) 


a 


r;  (2n  •*•  1)  \y^  Ya^kr)  *  c'a  ^^(kr)']  P^(cos  6), 


in  which  h  and  c  have  to  he  determined  from  the  boundary  conditions, 
n     n 

This  form  is  taken  since  the  function  u  results  from  the  superposition 

of  waves  reflected  from  the  upper  houndaiy  and  those  reflected  from  the 

lower  "boundary.  The  first  of  these,  which  are  incoming  waves,  are  seen 

from  (23a)  to  correspond  to  terms  of  the  form  h  "Xy  (kr).  The  waves 


^Stratton,  J.A.  :  Electromagnetic  Theory,  p  UU 


11. 


reflected  from  the  lov/er  toxindaxy  are  out-going  v;aves,   and  hence,   from 
(23'b),    correspond  to   terms  of  the  form  c    ^  (kr). 

Returning  to  the  case  in  which  the  refractive  indes  varies  in- 
versely as   the  square  of  r  in  layer  I,   v;e  can  also  vnrit©  the  Hertzian 
function  in  this  layer  in  the  form 

3y  applying  the  inversion  to  equations   (23a),    (23"b)   and  (2U),   ue  find 

^00  2  ,2 

^  =  I    ^  (2n  *  1)  ^i/^i^-)   ^J^)     P^Ccos  e).        r  >  '0.        (26a) 

.00  .2,2 

=  1    J^  (2n  +  1)  ^]/J^-)  ^^(:^)  P^Ccos  6).  r  <  '0.        (2Sb) 


v2 


u 
a 


I    ^2n  *  1)     [l>,\l/,(^)  -  c^  ^^(^)]    Pjcos  e)  .  (27) 


In  the  outermost  region,    since  the  refractive  indez  is  constant 
and  the  waves  are  out-going  only,   the  Hertzian  function  is 

00 

^n  "     Z  ^2a  •«•  1)  a     ^„(1:  r)     P  (cos  G)   .         r  >  d.  (2S) 

0         *-^  n    ^n     on 

In  the  earth,   also,   the  refraxstive   indez  is   constant,   and  the  waves  are 
in-coming  only.     Hence  we  have  for  the  Hertzian  fanction  inside  the 
earth, 

CO 

^i  =     H  ^2n  -i-  1)  d^  ^a(V^     P^(cos  G)   .  r  <  a.  (29) 

We  now  apply  the  houndary  conditions  at  r  =  d.  Since  £  u  must 

"be  continuous,  the  suja  of  (26a)  and  (27)  must  he  equal  to  (2S)  when  r  is 

2  k 
replaced  hy  d.  Writing  this  equality  and  dividing  "by  k  a  ,  we  have 


k2 


°   -n^n^V)=  I  l\^n(^)  ^^n  Xn^^^ 


k2 


"^Yn  ^^^  XJ'^   )1        (30) 


12. 

I^^        ^  ?  (•'-a) 

where  K  =  — ^  .     Siasilarly,  cinca    — i^6r~^    is  eontlimoiis,  we  Ealtiply 

(26a)  and  (27)  "by  r,   dlf f ereatiatie ,  replcco  r  "by  d  and  add,     Shla  result 
must  "be  equal  to  that  o'otaiacd  tj  leiltiplyijis  (2S)  "by  r,  differentiatiug, 
and  setting  r  =  d.     2hu3 

k  d     a     y  (^  ^)  »  -    -^(^    V(-^)  *  c     t'(^) 
o        nSn^o  ,2    La  Ta      an  *»a*  d     ' 

a 

In  obtaining  the  left-hand  neatsr  of  (3I)  ua  laaro  uced  tho  approsinatloa 
•4  t   ^n^V^l    r=d      =  V  Xn<^o^^    "^   ?n^-o^> 

since  the  tera  neglected  is  ssall  in  coiaparisoa  tjith  tho  one  retained. 

We  apply  the  "boundary  conditions  at  r  =  a  in  a  lil:s  Eannor, 

Since   £  u  is  continuous,  the  sua  of  (2613)  and  (27)  Eust  "be  e^»ial  to  (^) 

2 
whan  r  is  replaced  "by  a.     Writicg  this  equali^  and  dividing  "by  I:  , 

VQ  have 


2 

Since    ^^    ■    is  continuois,  we  nultiply  (sST))  and  (27)  Ijy  r,  differentiate, 

set  r  =  a,  and  add.     The  result  laust  "be  equal  to  that  ol)taiasd  ty  differ- 
entiating (29)  and  setting  r  =  a.     Thus  .  /. 

^^n  X'rSh-^  =  -  ^^  \.\  il^^^  *  =n  ^n^^->    "^  Yn^T^^   ^'n^^^     •       ^33) 

Equations(3o)  -  (33)  fona  a  systeia  of  four,  equations,  linear  in 
a  ,h  ,c  and  d  .  Solving  for  these  quantities  "by  detenainanto. 


13. 


'  we  olstaia 


a  = 


'^       ^ 


~n     a  Sn     a 


Y.W)   t;<!-)    --^M-)  ^^c^a) 


■a^'-? 


Vie 
ak. 


\|/,(V) 


"bk 


I  K'V' 


vi 


-V.(%-)^(S|.)       ^(£|., 


-Y.<T^)?r;'^'       ^>) 


-dk'= 


.T^k  d) 


br 


k  d        , 
-|~^(kd) 


l)k2 


YA-^ 


ak,  / 


A 


Ya<^)         -¥n(^)  ^.(^) 


dk 


•//ka 


y;(ka) 


ka  \  c-  /ke. 


-V.'^)4<^) 


/ka.' 


-¥.'^)  K.^^^ 


-■^J^^K>-^ 


■  bk' 

i4n 

0 

' 

k  d 
o 

•b£ 

[k^d) 

0 

-ak^^ 
bk^ 

0 

ak. 

bk 


\|/,(k,a) 


l|/^(k,a) 


lU. 


vk 


In  a    ^n  a 


1)Z 


k  d 


T-   rJV)   ^  -1^)  t(^) 


.^2    ^n^  0 


o""    /'  (k  d) 


n^  d 


n  0 


in  whicli 


A= 


^  n  d    ^n  a  ' 


Yn<^^   <^n^^^^ 


Y^Cka)   ^;(ka) 


-d  k 


■bK 
k  d 

0 


^ 


S  2.   O 

;-'(k  d) 

^  n^  0 


0 
0 


-ak 


3_  Y,  (v> 


^' 


^  if;(v) 


tk 


If  these  values  of  "b  and  c  are  inserted  in  (26a),  or  (2D'b),  and  27, 
and  th.e  resulting  expressions  for  u.  and  u  are  added,  we  have  the 
Hertzian  potential  function  at  any  point  in  the  layer  of  variable  index. 
Similarly  inserting  the  values  of  a  in  (2S),  or  d  in  (29),  we  have  the 
expression  of  the  Hertzian  function  at  any  point  in  the  outer  region,  or 
any  point  inside  the  earth.  It  is  not  hard  to  see  that  these  series  do 
not  converge  rapidly.  For  large  values  of  the  argument  the  functions  r  n. 
ajid  y  are  closely  approximated  hy  their  asymptotic  expressions,  and 
these  values  change  slowly  with  n.  Hence  succeeding  values,  "both  of  the 
coefficients  and  of  the  variaole  factors  in  the  series,  change  slowly. 


15. 


It  vould  therefore  "be  necassary  to  use  a  large  mim"ber  of  terms  to  ©"btaia 
a  good  approximation,  to  u.  The  computation  of  the  coefficients  of  these 
terms  would  clearly  entail  a  vast  aaount  of  labor,  which  does  not  seem 
justified  in  vieu  of  the  limited  application  of  the  results. 
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